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The influence of a small perturbation on a cavity mode plays an important role in fields like optical 
sensing, cavity quantum electrodynamics and cavity optomechanics. Typically, the resulting cavity 
frequency shift directly relates to the polarizability of the perturbation. Here we demonstrate that 
particles perturbing a radiating cavity can induce strong frequency shifts that are opposite to, and 
even exceed, the effects based on the particles’ polarizability. A full electrodynamic theory reveals 
that these anomalous results rely on a non-trivial phase relation between cavity and nanoparticle 
radiation, allowing back-action via the radiation continuum. In addition, an intuitive model based 
on coupled mode theory is presented that relates the phenomenon to retardation. Because of the 
ubiquity of dissipation, we expect these findings to benefit the understanding and engineering of a 
wide class of systems. 

PACS numbers: 42.82.Fv, 78.67.-n, 42.25.Hz, 42.25.Fx 


The fact that a small perturbation of a potential can 
influence the distribution of a systems energy levels is 
a well-known principle permeating various branches of 
physics. In quantum mechanics, for example, the effect 
of a perturbing potential H' on an eigenstate |' 0 o) is that 
it modifies its unperturbed energy Uq by an amount 5 U = 
where \ fj) is the new eigenstate and we assume 
(' 0 o|' 0 ) ~ 1 . In electrodynamics, a local change of poten¬ 
tial (z.e. permittivity) can impact the frequency of a 
resonant cavity. This is at the basis of many applications 
that use the influence of a perturbing atom, molecule, or 
dielectric body to establish an interaction that can be ex¬ 
ploited for optical sensing or control [THH] . The shift of a 
mode’s complex eigenfrequency cj = cjc — i/^/2, with cav¬ 
ity resonance frequency uJc and linewidth n, due to a local 
permittivity perturbation Ae contained in a volume AV 
is given hy 600/00 = — f^y dV [eoAeEo • Ep] / 4 Uo. Here 
Eq and Uq represent the field and total energy of the un¬ 
perturbed cavity mode, respectively, and Ep is the per¬ 
turbed field [illO]. In particular, if AV is small enough 
such that the perturbing particle can be described in 
the dipole approximation, the complex frequency shift 
is directly related to the particle polarizability a, read¬ 
ing 600/00 = —(a|Eop/ 4 [/o [ 9 ]. However, this widely em¬ 
ployed result of Bet he-Schwinger perturbation theory m 
is strictly only valid when radiation to the far field is 
negligible. In several recent developments, radiation loss 
proved decisive in determining a system’s eigenmode, e.g. 
for so-called ‘states bound in the continuum’ [IIHISl and 
in describing complex plasmonic resonators [m [ 15 ]. In 
that context, the question arises to what extent the con¬ 
ventional paradigm to determine perturbed cavity fre¬ 
quencies holds in practical, open, systems. 

Here, we study the eigenfrequencies of a radiating 
optical cavity as it is perturbed by carefully designed 
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resonant plasmonic nanoparticles. For a resonant per¬ 
turbation with center frequency 00g, and linewidth 7, 
the value of a strongly depends on the detuning A = 
<^c — ^a- In absence of radiation, the cavity mode thus 
redshifts(blueshifts) at negative (positive) values of A and 
the linewidth broadens near A = 0, in direct response to 
Re[a] and Im[(a], respectively (Fig [^). Tuning oog, by 
varying the length of the (plasmoni^ resonators allows 
to systematically study the induced cavity response as 
a function of detuning. Importantly, our experiment is 
designed such that radiation losses from the cavity mode 
and plasmonic resonators overlap. This, as we will show, 
leads to a strong additional contribution to 600 which is 
not captured by a. We study this ‘radiation interaction’, 
as we will call the effect, in detail and show that it can 
induce strong eigenfrequency shifts that are opposite to, 
and even exceed, the effects based on the particles’ po¬ 
larizability. 

A sketch of the system is shown in Fig.[^. The exper¬ 
iments are performed using a fundamental cavity mode 
(194.4 THz, Q^6.5xI0^, TE polarized) of a toroidal sil¬ 
ica microcavity fabricated on the edge of a silicon chip 
(~36 fim diameter, Eig. §) [isiin]. The cavity is per¬ 
turbed by gold nanoantennas deposited on a glass sub¬ 
strate (Eig. IT)’ which are controllably placed in the 
evanescent field of the cavity. The antennas (nanorods of 
length L, width 120 nm and thickness 40 nm) are aligned 
with their (long) principal dipole axis to the polarization 
of the cavity mode. A frequency-swept narrowband laser 
source (^ 0.7 fiW) is coupled into a tapered fiber that 
is brought close to the cavity. The transmission spec¬ 
trum through the fiber shows a Lorentzian dip around 
the cavity resonance frequency (Eig. [^), from which we 
determine the (perturbed) resonance frequency oOc and 
linewidth k of the cavity mode. Independent normal- 
incidence transmission measurements (Eig. [^) on the 
fabricated antenna arrays, in absence of the cavity, yield 
the normalized cavity-antenna detuning A/7. Notably, 
the presence of the glass substrate allows the cavity mode 
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FIG. 1. a) Cavity frequency and linewidth change due to 
perturbation by a small resonator, considering only the res¬ 
onator polarizability. Dashed lines indicate the real (top) and 
imaginary (bottom) part of the polarizability as a function of 
A/ 7 , b) An array of perturbing gold nanoantennas is placed 
in the near field of a toroidal microcavity. Light transmission 
through the tapered fiber is detected to determine the cavity 
eigenfrequency shift, c) Scanning electron microscopy (SEM) 
image of a high-Q silica microtoroid at the edge of a chip. 

d) SEM image of gold nanoantennas on a glass substrate. 

e) Transmittance through the tapered fiber around a cavity 
mode resonance. A Lorentzian fit (black solid line) is used 
to determine the optical linewidth (?^30 MHz) and frequency 
of the cavity. The sidebands result from a 150-MHz phase 
modulation used to calibrate the frequency axis, f) Normal- 
incidence transmission spectra of gold nanoantenna (design 
length L) arrays. Lorentzian fits (sold lines) give the antenna 
resonance frequency and linewidth. The dashed line indicates 
the cavity frequency. 


to radiate into the glass at a well defined angle just be¬ 
yond the critical angle and, similarly, allows scattering 
of antenna radiation into the substrate [18]. As the an¬ 
tennas are coherently excited by the cavity field, from 
which also the cavity radiation originates, the zeroth- 
order diffraction by the array is expected to overlap with 
the cavity radiation. 

To measure the change of the complex cavity reso¬ 
nance frequency due to the nanoantennas, we compare 
the frequency and linewidth of the cavity mode with and 
without the antenna array placed in the near-held of the 
cavity. Our measurement procedure m ensures that in 
both cases the sample is positioned at the same distance 
to the cavity. Figure shows the results for two differ¬ 
ent arrays of nanoantennas: In the top hgure, an antenna 
array {L = 360, pitch 800(900) nm along the long(short) 
axes of the antennas) is seen to induce a broadening of 
the cavity linewidth, together with a blueshift of the fre- 



FIG. 2 . a) Depending on A, antenna arrays induce linewidth 
broadening (top) and narrowing (bottom), while simultane¬ 
ously inducing the cavity resonance frequency to blueshift. 
The blue and red horizontal bars represent the linewidths of 
htted Lorentzian lineshapes (black lines), b) Shift in reso¬ 
nance frequency and linewidth of the cavity mode, normalized 
to its linewidth in absence of the antennas [19] , due to pertur¬ 
bation by antenna arrays with a constant pitch size and vary¬ 
ing antenna length. Top: For all detunings a blueshift of cav¬ 
ity resonance frequency is observed. Bottom: Approaching 
A = 0 induces linewidth narrowing. Both the blueshift and 
linewidth narrowing are contrary to the expectation based on 
the particles’ polarizability (Fig. [m). Error bars depict stan¬ 
dard deviation, vertical error bars fall within the plot markers. 
The solid lines represent a coupled mode theory fit. 


quency. The latter is clearly surprising, given the fact 
that A ^ — 1.27 for this sample m, where one expects 
a redshift if one only considers the resonant particle’s 
polarizability (Fig. &)■ The bottom panel in Fig. 
shows that for another array {L = 400, pitch 800(1100) 
nm) with slightly positive detuning (A « O. 87 ) again a 
blueshift is observed, but this time accompanied by nar¬ 
rowing of the cavity mode. Such a reduction of damping 
can obviously not be ascribed to the particle’s polariz¬ 
ability alone, as its imaginary part is necessarily positive 

(Fig#). 

To systematically study these effects, we gradually 
tune the resonance frequency uJa by varying the length of 
the antennas (Fig. HD- The array pitch is kept fixed at 
800(1500) nm along the long(short) antenna axes, chosen 
such that inter-antenna coupling is negligible. Moreover, 
in all examples we show, Bragg-scattering between clock¬ 
wise and counter-clockwise modes [ 20 ] is small enough 
such that mode splitting induced by the nanoantennas is 
smaller than the cavity linewidth. Figure shows the 
resulting cavity frequency and linewidth changes, nor¬ 
malized to the cavity linewidth on glass to allow averag¬ 
ing multiple scans m, as a function of A/ 7 . For these 
negative detunings, a consistent increase of the resonance 
frequency (i.e. blueshift) is observed that slightly rises 
as the antenna frequency approaches that of the cavity. 
As such, both the sign and the trend of 6 uJc are incon¬ 
gruous with the expectation based on the polarizability 
(Fig. ^). Moreover, the broadening of the cavity, ob- 
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served for a detuning of A — O. 67 , quickly reduces as 
the - strongly scattering - antennas are tuned closer to 
the cavity resonance. We note that these trends are con¬ 
sistently observed also at other periodicities [19] and that 
the shifts cannot be explained by thermal heating, neg¬ 
ligible at the employed powers. 

To explore the origin of these surprising results, it is 
necessary to consider the complete Bet he-Schwinger cav¬ 
ity perturbation equation m, derived without neglecting 
radiation: 


Suj = —uj 


a|Eo|2 

4C/o 


4C/o 


/ ' 


dA [((5E X H*) • h + (E* X 6 H) • n]. 


( 1 ) 


This full expression, which is also valid in any open, non- 
Hermitian system m, contains an additional integral 
term involving radiated fields at a surface 6 V (with nor¬ 
mal unit vector h), enclosing the same volume that was 
used to evaluate the total energy Uq. Here, SFi and 61i 
are the difference between the fields in the presence (E^, 
Hp) and absence (Eq, Hq) of the perturbation, and can as 
such be associated with the field scattered by the pertur¬ 
bation. Thus, the cavity eigenfrequency is additionally 
modified by an energy fiux that is evaluated by combining 
the fields of the perturbed and unperturbed eigenmodes; 
z.e, an overlap of scattering by the perturbation and di¬ 
rect cavity radiation. Remarkably, the integral is inde¬ 
pendent of the distance at which 5V is chosen, consistent 
with the fact that it can be associated with radiation of 
the system. To date, this far-held contribution has been 
omitted in practically all analyses of cavity perturbation. 
In a select number of experiments, a reduction of cavity 
linewidth was observed [2TII23] and tentatively attributed 
to interference of radiation of the cavity and scattering 
by the perturbation. Importantly, the integral suggests 
that in principle also the resonance frequency (i.e., the 
real part of 5uj) can be affected by the same mechanism, if 
the phase difference between cavity and scattered radia¬ 
tion in the far held (captured in Eq and (5E, respectively) 
dihers from 0 or tt. So far, the frequency has been ex¬ 
pected to rely only on the local variation of the applied 
potential jS] UHl [23H26] as it is contained in the hrst term 
of 0 . 

To verify the importance of radiation in a complete 
description of the induced cavity response, we perform 
hnite-element method eigenfrequency studies with and 
without antenna m in a similar two-dimensional sys¬ 
tem (Fig. |§), from which we extract the helds Eq, Hq, 
^E and ^H. This allows the calculation of the individ¬ 
ual terms of displayed as a function of detuning in 
Fig.|^. The hrst term (black crosses, ‘A’), contributes to 
Suj according to the expected behaviour for a perturbing 
resonator sketched in Fig. [TJi: it produces a dispersive 
detuning dependence of SuOc and a dissipative trend for 
Shz. This immediately disqualihes this term as an expla¬ 
nation for the measured changes in cavity resonance and 
linewidth. The second term, i.e. the contribution to the 



FIG. 3. a) Electric held prohle in a simulated system contain¬ 
ing cavity, antenna and substrate. Both antenna and cavity 
radiate into the glass under approximately the same angle, 
resulting in a radiation interaction. The contribution of this 
interaction is calculated on the dashed black line (integral of 
§). The arrow points to the position of the antenna, b) 
Top(bottom): the change in resonance frequency (linewidth) 
calculated using 0 . Black crosses (‘A’) show the contribution 
of the polarizability to the eigenfrequency shift Suj, which re¬ 
sembles the lineshape as we expect it from a resonant particle 
(Fig. [^). The radiative contribution (green stars, ‘B’) shows 
distinct blueshifts and linewidth narrowing near A /7 = 0. 
Adding both terms (blue circles, ‘A+B’) yields the complete 
cavity eigenfrequency shift, which matches the experimental 
trends. 


change in eigenfrequency related to the radiation inter¬ 
action (green stars, ‘B’), shows a dramatically different 
behavior: it causes distinct blueshifts and linewidth nar¬ 
rowing around A /7 = 0. Surprisingly, the magnitude 
of the blueshifts induced by this term can even exceed 
the contribution due to a. This is a key point of our 
observation, stressing the importance of this new contri¬ 
bution in a model system. The sum of both terms (blue 
circles, A+B) yields the complete eigenfrequency shift of 
the perturbed cavity, which now qualitatively matches 
our experimental results; blueshifts and linewidth nar¬ 
rowing around A /7 = 0. These calculations thus con¬ 
firm that radiation interactions in open systems can lead 
to both cavity blueshifts and linewidth narrowing and 
importantly, that this effect can even dominate the total 
change in eigenfrequency of a cavity mode. 

It is essential to realize that a modification of the cavity 
resonance frequency (in contrast to the alteration of the 
linewidth) due to back-action via the radiation contin¬ 
uum only occurs when the phase Scf) between cavity and 
antenna radiation is different from 0 or tt (Fig. [^). The 
exact strength of this back-action, and as such the sum of 
direct {i.e. due to a) and indirect {i.e. via the radiation 
continuum) coupling between cavity and antenna, will 
depend on the overlap between their radiation profiles 
and on the value of Scj). 

This phase difference and its physical connection to 
back-action can be captured in a simple model (Fig. [^) 
based on temporal coupled-mode theory m- Let us as¬ 
sume for now that the cavity and antenna are coupled 
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FIG. 4. a) The cavity and antenna are coupled with rate pc, 
which is the near-held interaction scaling with the polariz¬ 
ability of the antenna, and can radiate into the continuum 
with a certain delay ((50) with respect to each other. The 
interference between both radiation prohles and the overlap 
with the cavity mode results in a back-action effect (single¬ 
headed dashed arrow) on the cavity mode, b) Coupled-mode 
theory. The cavity and antenna are coupled with rate pc and 
couple to the continuum with complex rates = p^ and 
9i — 9i respectively. The difference between (j)^ and 0^, 
summed with the phase response of the antenna, is now equal 
to the phase difference (50. /^o and 70 are the losses of cavity 
and antenna into non-coupled channels. The environment is 
modeled as an oscillator with decay rate 


with complex coupling rate pc = pc and can both in¬ 
teract with the environment (glass) with coupling rates 
9 k = 9 i<i and p^ = p^ , respectively. We treat the 
joint environment as a separate mode to allow backac- 
tion, taking the limit of large decay rate 0 to mimic a 
broadband continuum. The cavity and antenna, consid¬ 
ered separately, each decay via this continuum at a rate 
Hi = 4|p/^p/0 and 71 = 4|p^p/0, respectively. This is 
distinguished from decay into all other, non-overlapping, 
modes at rates hq and 70 , such that hq ^ Hi = h and 
7 o + 7 i = 7 - Importantly, the phase difference between 
cavity radiation (blue path. Fig. Ef) and antenna scat¬ 
tering (orange path. Fig. [^) in this model now directly 
relates to (50 via (50 = tt + — arg[A + 17 / 2 ], where 

= 0c+ 0/^ + 0^, and arg[A + i 7 / 2 ] is the phase response 
of the antennas m- The natural mode of interest, with 
complex frequency cj, of this hybrid system can be found 
by equating the determinant of the coupling matrix [19] 


/oj-oJc + iKo/2 g* gi^\ 

M = \ gc A + i7o/2 g* 

\ 9i 57 '<n) 

to zero. Solving for complex uj yields 

5u} - \5 k/2 = ~ ^i7i/4 + igcV/ti7i cos^ 


( 2 ) 


(3) 


which depends on the radiation overlap /^i 7 i, (constant) 
phase and the coupling rate Pc- Because the over¬ 
lap and phase cannot be independently chosen, it is 


impossible to directly fit © to the experimental data. 
Therefore, we retrieve (50 from the finite element simu¬ 
lation and, using the relation between 5(j) and ff>, obtain 
= 3.76 X 10 “^ 7 r. Note that this value of implies 
that ( 50 ^ 7 r/ 2 atA = O, completely opposite to the con¬ 
ventionally assumed case where destructive interference 
(J 0 = tt) only contributes to the linewidth of the cavity 
mode. Fixing at this value and constraining the ex¬ 
perimentally found values of /i:, A and 7 (retrieved from 
independent spectroscopic measurements. Fig. [^ , a fit of 

g yields excellent correspondence to our data in Fig. If 
•lid lines), resulting in an overlap of Hi^ijn^ = 0.68 
and coupling rate Pc/27r = 21.6 GHz. 

Concluding, we have shown how coupling through 
the radiation continuum leads to cavity blueshifts and 
linewidth narrowing in a coupled cavity-antenna system. 
Surprisingly, the non-local, radiative effect on the cavity 
is larger than the induced cavity response due to local 
perturbations, and is not just an artefact of large radia¬ 
tion losses but can even manifest itself in a high-Q cavity 
such as studied here. Therefore, similar effects are ex¬ 
pected in many other systems where cavity and perturba¬ 
tion can radiate into the same channels, such as photonic 
crystals |2TH23l|28] and Fabry-Perot etalons [2511^ . In 
addition, it was shown that there is a direct link between 
the phase (50 in the Bet he-Schwinger equation and the 
contributions to back-action in a coupled-mode model. 

Interestingly, these results could also shed new light 
on radiative interactions in strongly radiating systems 
such as metamaterials and complex plasmonic res¬ 
onators. In El, it was postulated that radiation can 
be included in the coupling between two resonators 
by using a complex coupling rate, while in [32|, a 
complex coupling rate was argued to originate from a 
complex-valued extinction cross-section. In fact, such 
a complex coupling rate between two resonators is 
easily derived from the coupled-mode model presented 
in this work [19]. Furthermore, it will be interesting 
to see how the concepts and theory established here 
relate to recently developed methods to relying on the 
normalization of leaky modes to describe the response of 
complex photonic systems [33][34]. In a different context, 
recent advances in optomechanics using simultaneous 
dispersive and dissipative (z.e radiative) coupling [35ll37] 
have been reported to enable on-resonance cooling [35] 
and sensitive readout of nanomechanical motion [36] . 
As such, we expect our work to be directly important to 
the design of novel optical devices such as sensors, for 
example, and benefit the understanding of the physics 
of open systems in optics and beyond. 
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Foundation for Fundamental Research on Matter (FOM), 
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from an NWO Vidi grant. 
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Supplementary Material 


I. BETHE-SCHWINGER CAVITY PERTURBATION IN OPEN SYSTEMS 


In this section, we show the derivation of the exact Bet he-Schwinger cavity perturbation formula, analogous to the 
approach by mm- Afterwards, we discuss how small approximations result in equation 0 of the main text. We note, 
that the calculation of the different contributions to the eigenfrequency c hang e in our simulations did not involve the 
approximated equation Q, but the exact formula as shown in equation (S19). 


A. Eigenmodes 

We consider the modes of an open cavity described by a spatial distribution of permittivity eoe(r) and permeability 
jiiojii{r) (in the following, the spatial dependence of both is implicitly assumed). The eigenmodes of the system are 
found by solving Maxwell’s equations in all of space in the absence of external drives. This yields solutions of the 
form 


E(r,t) = Eo(r)e-‘“*, (S4) 

H(r,t) = Ho(r)e-‘‘"*, (S5) 

which satisfy Maxwell’s equations for some complex frequency 

cj = cjc-i^, (S6) 

where ujc and hz denote the (real) resonance frequency and energy decay rate, respectively. Typically, we will consider 
systems where hz uJc. Important relationships between Eq and Hq include 

V X Eo = icj/io/^HoV X Eg = —icj/io/^Hg, (S7) 

V X Ho = —icjeoeEoV x Hq = icjeoeEg. (S8) 


B. Dielectric perturbation 

We now consider that the cavity is perturbed, such that in a finite volume AV, the relative permittivity and 


permeability are changed to new values 

ep = e + Ae, (S9) 

= /^ + A/i. (SIO) 

This perturbed system has different eigenmodes written as 

E(r,0 = Ep(r)e-‘“>’* = (Eo(r) + 5E(r)) (Sll) 

H(r,0 = = (Ho(r) + <5H(r)) (S12) 

where Ep(r) and ujp describe the spatial dependence and resonance frequency of the perturbed cavity mode. Defining 
SEi = Ep — Eo, SH = Up — Ho and 5uj = ujp — uj yields 

V X (Eo + SIEj) = i (cj “h Siu) f^ofip (Ho + ^H), (S13) 

V X (Ho (^H) = —i {uj -|- 5uj^ (Eo <^E). (S14) 

Combining the above equations for the curl of the fields gives 

V X ^E = icjp/io/ipHp — icj/io/iHo, (S15) 

V X (5H = —iujpeoep'Ep + icjeoeEo. (S16) 


We next take the dot product of Hq and Eq with the curls of SEi and ^H, respectively and rewrite both making use 
of the vector identity 


a • (V X b) = b • (V X a) — V • (a X b). 


(S17) 
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Subtracting the obtained expressions gives us 

6lu (cocEq • Ep + /io/^HQ • Hp) = — (cj + 6uj) (cqAcEq • E^ + /ioA/iHg • H^) 

- i (V • (^E X HS) + V • (ES X 6 li)). (S18) 

Finally, we take the integral of both sides over a (very large) volume and apply Gauss’s theorem to arrive at the 
Bethe-Schwinger equation [3]: 

Ju; [ dV [eoeES • E^ + /io/iHS • H^] = - (u; + (^cj) [ dV [eoAcE^ • E^ + /ioA/iH^ • H^] 

Jv Jav 

-if dA[(^ExHS)-h+(ES X(5H)-h], (S19) 

Jdv 

where dV denotes the boundary of V and h the unit vector normal to that surface. 


C. Approximating the cavity perturbation 


The above equation (S19) yields an exact expression for the (complex) cavity shift 

S tv 

Scj = Scjn — i—. 

2 


(S20) 


To arrive at equation 0 of the main paper, we start by making the following approximations: If the effect of the 
perturbation is small, then cj-hSu; ^ cj, and when the total volume is large compared to the volum e of t he perturbation, 
it will be allowed to omit the contribution of Eg • dFi in the integral on the left side of equation (S19), which can now 
be written as 


6 uj 


[ dy[eoeES-E^ + /io/iHS-H^] 
Jv 


Suj [ dV 

Jv 


eoe |Eo|^ +/uo/i |Ho|^ = ASuUq, 


where Uq denotes the energy stored in the unperturbed cavity. Now the Bethe-Schwinger equation becomes 

6 lo = -^ [ dV [eoAeES ■ Ep + ■ HJ - ^ / dA [(,5E x HJ) ■ n + (E^ x SB.) ■ n]. (S21) 

4^0 Jav 4^0 Jav 

If the perturbation is caused by a variation of the permittivity (A/i = 0) and if the dimensions of its volume AI^ are 
small compared to the scale over which the cavity field Eg varies, one can make a useful approximation by considering 
the perturbation as a polarizable (point) dipole. Its total dipole moment p is given by 

where we introduced the dipole moment per unit volume, which we identify as 

^(r) = eo {cp - e) Ep(r) = eoAeEp(r). (S23) 

Taking into account that Eq can be considered constant over AV, we can now write the integral over AV in equation 
as 


ES • eo / . 

Jav 


d^AeEp = ES • p = ES-Of Eg, 


where we have introduced the polarizability a defined in this case through 

p =(a Eg. 


(524) 

(525) 


If a is diagonal, and we denote the magnitude of the polarizability in the direction of the cavity field by the scalar a, 
the integral over AI^ reduces to a |Eg|^. In this approach, we arrive at a new version of the Bethe-Schwinger equation 
for the perturbation of the cavity frequency: 

Suj =- ^ [ dA[{SEx-tia-n+(E*ox5U)-n], (S26) 

4(7o 4(7o Jav 

Here we see how the real part of the particle’s polarizability contributes to a frequency shift of the cavity, but so 
does the imaginary part of the surface integral. Likewise, the change of the cavity linewidth is related to the imaginary 
part of a and the real part of the integral over dV. 








Detector 



FIG. S5. Experimental setup. See text for details. 


II. EXPERIMENTAL DETAILS 

A. Sample fabrication 


Gold nanoantennas are fabricated in an array on a 170 /im thick glass coverslide. To start, a 130 nm layer 
of ZEP-520 resist is spin-coated on top of the coverslide. The nanoantennas are patterned into the resist using 
electron-beam lithography. After development, thermal evaporation of gold and a lift-off step yield the desired 
antennas. The antenna width and thickness were designed to be 120 and 40 nm and the length was varied between 
400 and 360 nm. The pitch along the short-axes of the antennas was varied between 750 and 1500 nm, and the pitch 
along the long-axes was fixed at 800 nm. A high-Q silica microtoroid (diameter « 36 fim) is fabricated on the edge 
of a silicon sample, largely following methods as previously reported In this work, spin-coating (ma-N 2410) 

and subsequent cleaving of the sample enabled targeted e-beam lithography of the disks. 


B. Experimental Setup 


Figure shows the experimental setup used in this work. We use a tunable fiber-coupled external cavity diode 
laser (New Focus TLB-6728, <100 kHz linewidth) to probe the microcavity. The laser frequency is controlled by 
an external voltage from a function generator. Coupling of light into the cavity is achieved using a tapered optical 
fiber, of which the position is controlled using piezo positioners (not shown). The electro-optical modulator (EOM, 
EOspace) is used to generate sidebands of the cavity mode at known RE frequency, which allows for calibration of 
our optical frequency axis. Eiber polarization controllers (EPC) ensure effective coupling into the optical cavity mode 
of interest. Before detection, the light that is transmitted through the tapered fiber is amplified using an Erbium 
Doped Eiber Amplifier (EDEA). Radiation leakage from the microcavity is collected using an objective (Nikon, CEI 
Apo TIRE lOOx) to obtain information regarding polarization and mode profile (through fourier-space imaging and 
polarization analysis of the back focal plane (BEP), Eig. [S^ of the cavity mode. Using this method we identify a 
fundamental cavity mode, of which the polarization is aligned with the long axes of the antennas. Due to limitations 
with regards to the numerical aperture (NA) of our objective (NA 1.33), imaging of the complete cavity mode-profile 
is impossible, thus resulting in a cut-off of the imaged mode profile. Considering this cut-off at NA=1.33 and the 
refractive index of silica (1.46), we estimate an effective mode index of around 1.35. We check the frequency stability 
of the laser itself using a fiber loop reference cavity (ELC). 
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FIG. S6. Mode characterization. Illumination of the BFP of our objective with white light (not shown) is used to calibrate 
the image on our infrared camera. This involves fitting the NA = 1 ring and maximum collection angle of our objective (orange 
circles), which are imaged on the camera. Sample tilt results in a slight offset between the outer and NA=1 ring. Subsequent 
monitoring of the radiation leakage of our cavity mode enables its characterization (fundamental TE mode) and an estimate 
(see text) of the effective mode index (?^ 1.35). The blue solid lines serve as a guide to the eye. 




walked path (steps) 

FIG. S7. Scanning along an array, a) Gartoon of a scan of the microcavity along the antenna array. Gareful alignment 
allows us to scan over the antennas by moving our microcavity between two predefined points (ref 1 and ref 2), which serve 
as our reference measurements on glass. Between those points we measure the effect of the antennas on the cavity mode, 
b) Evolution of the cavity linewidth during a scan along the antenna array, where each point is retrieved from a fit with a 
Lorentzian lineshape. When scanning along the antennas, we make sure that we see a plateau of data points. In our analysis, 
we only include the data points indicated by the grey regions, thus excluding the data taken during transit between reference 
and antennas. Eor the calculation of Sk, we average all the points clearly on the plateau and subtract the average reference 
value on glass (/^ref = l^snh + /^int). 


C. Measurement procedure 


Figure [S7p shows that the experimental data is obtained by stepwise scanning the microcavity, which is placed 
on a piezo-electric stage, along an antenna array and a stretch of bare substrate at a distance of ^1.1 /im (based 
on simulations). During each scan, the cavity-fiber distance is kept constant, which is checked before and after 
each experiment when the cavity is positioned at opposite positions adjacent to the antenna arrays. This check is 
performed by monitoring the cavity linewidth (broadened by the presence of the substrate), which depends strongly 
on the distance to the substrate. Each scan consists of 80-120 steps and multiple scans (2-5) are performed on each 
array. The change in resonance frequency {5uJc) and linewidth {5n) due to the antennas is obtained by subtracting 
the cavity resonance frequency and linewidth without the array placed in the near-held of the cavity, from the cavity 
resonance frequency and linewidth with the antenna array placed in the evanescent held of the cavity mode. Example 
data for the change in linewidth is shown in Fig. 873. To correct for any remaining sample tilt and/or slow drift, 
reference measurements on glass taken at the start and end of a scan are averaged. The values of 5uJc and are 
calculated by excluding data points taken during the transition from glass to array. 

Even though we check that all measurements are performed at the same cavity-sample distance, we would like to 
eliminate any errors due to small variations of that distance. To this end, we retrieve parameters that are robust 
towards small distance changes by considering the relative change in linewidth(resonance frequency) due to the 
antennas: 5ujc )/where /^sub is the broadening the cavity mode experiences due to the glass substrate (/^sub = 
z^ref ~ with the dccay rate of the bare cavity). Figure 88 shows that this method results in a relative change 
of the resonance frequency ujc and linewidth k which is independent of the distance to the substrate. It is obvious that 
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FIG. S8. Validation of dimensionless parameter. The measured effect on the cavity linewidth (a, bottom) and resonance 
frequency (b, bottom) due to the antennas is clearly increasing when the cavity travels a certain distance towards the antennas. 
As shown in the top panel of (a) and (b), our analysis results in a constant dimensionless quantity for the change in linewidth 
and frequency which does not depend on the distance to the antenna array. 


this quantity now allows straightforward averaging over multiple scans on the same array and comparison between 
different antenna arrays. Errors on the mean are calculated for 5uJc and while the error (one standard deviation) 
on uja (and thus the detuning) is retrieved from a Gaussian fit to the sum of squared residuals, obtained by displacing 
the fit result with respect to the normal-incidence tranmission spectra of the antenna arrays. 


III. FINITE-ELEMENT SIMULATIONS 


To gain more insight in our experimental results, and to calculate the different contributions to the complex 
frequency shift 5uj^ we perform 2D eigenfrequency simulations using the COMSOL software package (v5.0). The 
cavity (diameter = 36 fim) has a refractive index of 1.35, comparable to the experimentally determined effective index 
of our cavity mode and is placed 1.1 /im from the substrate, which has a refractive index of 1.5. The antenna is 
assigned a varying permittivity to mimic a change in detuning between cavity and antenna and for the surrounding 
we take a refractive index of 1. The simulation includes the use of a Perfectly Matched Layer (PML). 

In our experiment, changing the length of the antennas results in a changing plasmon resonance, in turn causing 
a varying detuning between cavity and antenna resonance frequency. To mimic this in our simulation, we assign a 
tunable permitivitty, which has a Lorentzian lineshape, to the perturbing antenna(s): 


S * 7s/2 


(S27) 


Sweeping the variable Ag now allows to calculate the change in cavity eigenfrequency as a function of detuning. The 
parameter S can be used to control the scattering strength of the antenna, while the antenna linewidth is held fixed 
at 7s = 1. To determine a sensible value of S, we analyze the eigenfrequency shift of the cavity when perturbed by 11 
antennas (pitch 1500, FigJS^) for different values of S and choose a value (S=340) which yields the best resemblance 
to our experiment (Fig. |S9|3). The values for the relative change in linewidth and resonance frequency are obtained 
via similar analysis as performed in the experimental situation, thus involving separate simulations for the bare cavity, 
and cavity perturbed by just a glass substrate. 

It is important to realize that the values Ag and 7 g are input parameters in the simulation, which do not include the 
effect of radiation losses on the antenna linewidth and resonance frequency [6] . To allow for a fair comparison with our 
experiment, we fit the results retrieved from the simulation with a Fano lineshape (black lines in Fig.|S9^), to obtain the 
radiation corrected values A and 7. Using these corrected values of A and 7, we rescale our horizontal axis and obtain 
Fig. [ 89 ) 3 . Note that this approach resembles the experimental situation, where the radiation corrected values of A 
and 7 are retrieved from the normal-incidence transmission spectra. It can be observed from Fig.|S9|3 that the orange 
(resonance frequency) and pink (linewidth) circles are a good, but not perfect, match to the experimentally obtained 
cavity frequency and linewidth change (blue circles). We attribute the discrepancy to the intrinsic difference between 
a 2D simulation, where the cavity is modeled as an infinite cylinder, and the experimental situation, which employs 
a toroidal cavity. This will necessarily change the radiation profile determining overlap and phase. Nonetheless, the 
2D simulation reproduces all the main features observed in the experiment. 
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-6 As/Ys ° A/y ° 

FIG. S9. Using multiple scatterers to find the scattering strength S. a) Filled circles: change in resonance frequency 
(top) and linewidth (bottom) of the cavity when perturbed by 11 antennas. The values on the horizontal axis are the input 
parameters in the simulation, which do not include effects due to radiation losses. Performing a fit (black line) using a 
fano-lineshape allows us to obtain the corrected values A and 7 , which are used in (b). b) The blue circles and solid lines 
are experimental data points and a fit using the coupled-mode model, respectively. The orange (top) and pink (bottom) filled 
circles are the relative frequency and linewidth change that are obtained from a simulation with S=340. The rescaled horizontal 
axis is obtained from the fit in (a). 


A. Cavity and perturbation without glass 

When calculating the complete Bet he-Schwinger formula, it is useful to first consider the situation where a cavity 
is perturbed by a single perturbing particle, in the absence of a substrate, as this should not allow for radiation 
interaction taking place, and thus result in a perturbation of the cavity eigenfrequency which is solely governed by the 
polarizability of the perturbing particle. The results of this simulation are shown in Fig. |SlQ^ . It is clear that the cavity 
indeed lacks radiation interaction (green stars) with the antenna and experiences linewidth and resonance frequency 
shifts that are solely related to the polarizability of the particle (black crosses), such that the total eigenfrequency 
change (open blue circles) is completely dominated by local effects. We remark that the Bet he-Schwinger cavity 
perturbation formula, which needs fields as it input, perfectly reproduces the cavity eigenfrequency shifts as predicted 
by the eigenvalues (orange(pink) filled circles) directly obtained from COMSOL. 


B. Cavity and perturbation with glass 

As discussed in the main text, introducing a glass substrate enables and extra coupling channel between antenna 
and cavity due to overlapping radiation profiles. Fig. |S1Q| shows results quantifying this effect for a cavity perturbed 
by a single antenna. Here the radiative interactions (green stars) are of importance when calculating the total (open 
blue circles) eigenfrequency shift of the cavity using the Bet he-Schwinger equation, dominating the effect of local 
perturbation (black crosses). Note that the calculation of complex Suj using the full Bethe Schwinger equation, 
yields almost exactly the same results when compared to the cavity eigenfrequency changes based on the eigenvalues 
calculated in the simulation (orange(frequency) and pink(linewidth) filled circles). We attribute the small difference 
that is observed mostly to the way joined data sets are handled in COMSOL. The use of joined datasets is necessary to 
obtain the fields ^E, but causes an error when integrating the closed-surface flux integrals (last terms of equation 0). 
The reason is that for joined datasets COMSOL projects full datasets on a less dense intermediate mesh [7]. We 
performed mesh convergence studies to see if this effect would diminish with increasing mesh density (up to 60 
elements/Ao). Although we observe that both the eigenvalues and the fields of our solutions are clearly converged, the 
small difference in the cavity eigenfrequency shift that we obtain using both methods does not disappear. Nonetheless, 
the discrepancy is small enough to reproduce all observed trends. 


C. Subtracting fields 

Figure]^ of the main paper shows the field profile Ep of a cavity perturbed by an antenna on a glass substrate. 
To increase the visibility of the radiation pattern, that figure and Fig. |S11| show the field amplitude for a cavity 
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FIG. SIO. Benchmarking the calculation, a) In the absence of a substrate, the cavity experiences a change in resonance 
frequency (top) and linewidth (bottom) that is solely associated with the polarizability (black crosses) of the perturbing particle, 
such that that the total eigenfrequency change as calculated by the Bet he-Schwinger equation (open blue circles) is completely 
dominated by local effects, leaving the radiative interaction (green stars) to be negligible. As expected, the Bethe-Schwinger 
cavity perturbation formula is able to perfectly reproduce the cavity eigenfrequency shifts when compared to the shifts obtained 
using the eigenvalues (orange(pink) filled circles) calculated by COMSOL b) In the situation where the substrate is present, as 
discussed in the main text, radiative interactions (green stars. A) are of importance when calculating the total (open blue circles, 
A+B) eigenfrequency shift of the cavity using the Bethe-Schwinger equation, dominating the local effect of the perturbation 
(black crosses, B). Note that the calculation of the full Bethe Schwinger equation yields almost perfect overlap with the cavity 
eigenfrequency changes as calculated by COMSOL (orange and pink filled circles), which are based on the eigenvalues of the 
solutions. We attribute the small difference that is observed mostly to the way joined data sets are handled in COMSOL, which 
involves projection on an intermediate mesh. This is known to cause loss of information and introduces a small error when 
integrating radiation patterns. 


placed 0.3 jam from the substrate and antenna 
find the unperturbed cavity field profile (Fig 


To calculate the full cavity perturbation formula, it is necessary to 

in order to obtain the field (iE 


_ „ _ . _ .Sll r, Eo) and subtract this from Ep 

(Fig. |Sll] 3). The teal line surrounding the field profiles represents a perfectly matched layer (PML), which is used in 
all calculations discussed in this paper, while the dashed black line is the boundary which we used for the calculation 
of the surface integral of equation 0 - 

A pitfall in subtracting fields from different eigenfrequency simulations lies in the arbitrary choice of amplitude and 
phase of the solutions. To overcome this problem, we monitor the amplitude and phase in both simulations (with and 
without antenna) in a point where modification of the field profile due to the perturbing antenna is expected to be 
negligible. This can be, for example, a point in the field maximum of the cavity mode far away from the perturbation, 
thus allowing to correct for possible differences in amplitude and phase between the calculations. To obtain the data 
in Fig [It, a single phase correction was performed, to match the phase at Ag = —3 to the unperturbed phase, in a 
mode maximum at large distance from the perturbing particle. 


IV. COUPLED MODE THEORY 

We will describe the interactions between cavity and antenna using coupled mode theory [8], essentially modeling 
the system as coupled harmonic oscillators with Lorentzian response. We note that while it is possible to treat the 
antenna as a harmonic oscillator without making the Lorentzian approximation, all essential physics is captured in 
the simpler Lorentzian model. We treat the decay in a common (overlapping) radiation continuum, and its possible 
backaction on both antenna and cavity, by considering the ‘environment’ as a third mode coupled to both antenna 
and cavity as well as to an independent decay port. The complex mode amplitudes of cavity (a), antenna (b) and 
environment (c) are normalized such that their magnitude squared equals the energy in each degree of freedom. They 
are collected in vector a = (a,6, c)^. The equations of motion of the system (in absence of driving fields) can be 
written as 


da 

dt 


(-iU-F)a, 


(S28) 
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FIG. Sll. Simulated field profiles. 

The field profile SEi of the perturbing 
Fig. of the main paper. In both (a) 
(PML). 


a) The simulated field profile Eo for a cavity in the absence of perturbing particles, b) 
antenna. This profile is obtained by subtracting the fields in (a) from those shown in 
and (b), the teal line surrounding the field profiles represents a perfectly matched layer 


where 


UJc 

-~9l 

-9^ 

-9c 


-~9: 

-~9l 

-h 

LUe 


(S29) 


is a Hermitian matrix containing the resonance frequencies and coupling rates. These coupling rates can in general 
be complex, and we will write them later as gj = with real gj and (pj. The matrix 

[^o /2 0 0 \ 

F = 0 7o/2 0 (S30) 

V 0 0 c/ 2 / 

is also Hermitian and contains the damping rates to the different independent coupling channels, which are associated 
with an individual mode a, 6, and c. As such, the output fields s = ( 5 ^, 55 , 50 )^ in these ports can be written as 
s = iFa, where 


0 0 \ 

K = [ 0 0 . (S31) 

V 0 0 vcj 


These equations satisfy time reversal symmetry and conservation of energy, assured by the hermiticity of Q and F 
together with the condition iF*iF = 2F. We now introduce the detunings Ac = cj — cjc, A = uJc — ujq and Ae = uJc—^e- 
The eigenmodes can be found by solving the system Ma = 0 where 


M = Iuj-VL-^iV 


/Ac + iKo/2 gt Qk 

I 9c Ac + A + i7o/2 

\ 9k 9i Ac + Ae + iC/2 


(S32) 


Since we assume the cavity and antenna are weakly coupled and we are only interested in the perturbed cavity mode, 
we can neglect Ac in the second and third row of M. The solution is 


Ac 


9I (A + + flc (^e + i|) - ‘^9c9k9i cos $ 


(S33) 


where we have introduced <F = 0c + + 07 - 

It is instructive to consider the case where the cavity mode is absent (a = 0^ gc = g^ = 0). In that case, the power 
radiated via the environment into port Sc is equal to (g‘^\b\‘^/{A^ + C^/4)- We therefore define the rate 

A2 + C74’ 


7i = 


(S34) 
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which is the rate at which the antenna decays into the common radiation continuum if the cavity were not there. 
Likewise, we consider the case without the antenna {b = 0^ gc = = 0), and define 


A2 + C74’ 


(S35) 


which is the rate at which the cavity would decay into this channel in absence of the antenna. In the latter case, the 
eigenmode would have detuning 


= (S36) 

where /i: = /^o + is the total decay rate of the cavity without the antenna. We now specifically write the general 
solution as Ac = Ac, such that 5Ac is the change of the eigenmode frequency due to the presence of the 

antenna and 5Ac = dujc — Taking the limit of C to infinity, which effectively amounts to setting Ae = 0, it reads 


(^Ac 


gl - + i^fcT^^/kycos^ 

A + i7/2 


(S37) 


where we have introduced g = (niji) / (nj); a number between 0 and 1 that is related to the overlap of the radiation 

of both modes. This is expression is equivalent to equation (§ in the main text. 

Next, let us turn to the relation of this model to the phase difference Sc/) as it appears in the Bet he-Schwinger 
equation. To this end, we express the field in the environment c, compared to the field in the cavity a. Its full 
expression (still taking Ae = 0) is 


c 2ie (A + if) - 

a - — -aTT|- 

The phase 5(j) must be equal to the phase difference of the above in absence of the particle {gc = ^7 = 0) and its phase 
in absence of direct coupling between cavity and environment {g^ = 0). Performing this subtraction yields 

(50 = TT + 4> - arg ^A + i^^ , (S39) 

showing that the phase difference of the two decay paths is directly related to the (relative) phases of the coupling 
rates via 4> and the phase response of the antenna [9]. It becomes clear that the standard case, where the two paths 
interfere destructively on resonance (A = 0), occurs for 4> = 7r/2, such that dcj) = tt. 

Finally, we note that it is straightforward to eliminate the environment mode c from the model to construct an 
equivalent system of equations M'{a^ b)^ = 0. In that case the total damping rates k and 7 of the uncoupled cavity 
and antenna appear on the diagonal: 


M' = 


■if 


(5c + “*’) 


e {gc + 57^1716’^) 


■ 

^2 


We can choose the phase (p^ = 0 without loss of generality, such that 


(S40) 


M' = 


gc + IfKilie 


-i<I> 


9c + jVKi7ie‘*\ 

A + q )■ 


(S41) 


Note that the off-diagonal matrix elements of M', representing coupling between cavity and antenna, are now no 
longer Hermitian. 


V. VARYING ANTENNA PITCH 


The text of the main paper shows the results and analysis of antenna arrays with pitch 800( 1500 ) nm along the 
long(short) axis of the antennas. More pitches are, however, experimentally investigated and Fig. SI2 shows the main 
results which we obtain when varying the pitch along the short axis (750-1500 nm) of the antennas, while the pitch 
along the antenna long axis is kept constant at 800 nm. In this figure, the points which are surrounded by a dashed 
circle correspond to the arrays used to generate Fig. of the main text. The main conclusion that can be drawn 
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FIG. S12. Varying the antenna array pitch. Changing the pitch along the antenna short axis (between 750 and 1500 nm) 
does not greatly affect the cavity eigenfrequency shift. For all arrays, we observe cavity blueshifts (top panel) and linewidth 
narrowing (bottom panel). Note that in the bottom panel, the different pitches cross the point of zero linewidth change at 
slightly different values of A/y. Modification of the single particle polarizability due to the array is the most likely cause of 
this behaviour. The values on the vertical axis are corrected for the unit-cell area and the points which are surrounded by a 
dashed circle correspond to the arrays used to generate Fig.[^ of the main text. Errors on the mean (vertical axis) fall mostly 
within the plot markers and the error on the horizontal axis is one standard deviation, retrieved from a Gaussian fit to the 
sum of squared residuals (obtained by displacing the fit result). Due to limited detector range, a reliable horizontal error for 
the data at P1100,A/7 0.8 could not be obtained using this method. The error on this data point was thus set equal to the 

error for array P1100,A/7 —0.15. 


from Fig. |S12| is that, qualitatively, the different arrays behave the same: all arrays induce blueshifts of the cavity 
resonance frequency (top panel) and linewidth narrowing (bottom panel) of the cavity on approaching A/y = 0. With 
changing pitch, the value of A/y where the ’’relative linewidth change” equals 0, slightly varies. We attribute this 
effect to the influence of the array on the single particle polarizability. Remarkably, for the largest possible detuning 
we measure linewidth narrowing up to 30% for a pitch of 1100 nm. 

As in the main paper, the values of A/y are retrieved from a fit of the normal-incidence transmission spectra using 
a Lorentzian lineshape. Lattice sum calculations uni are performed to exclude strong modification of A due to the 
change of effective angle of incidence (in the experiment) with respect to the normal-incidence transmission spectra. 
In addition, total internal reflection measurements on the samples with a 1500 nm pitch showed no significant change 
of A with respect to the values obtained from the normal-incidence measurements. 
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